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In polymer: fullerene solar cells, the origin of the losses in the field-dependent photocurrent is still controver- 
sially debated. We contribute to the ongoing discussion by performing photo-induced charge extraction mea- 
surements on poly(3-hexylthiophene-2,5-diyl):[6,6]-phenyl-C6i butyric acid methyl ester solar cells in order to 
investigate the processes ruling charge carrier decay. Calculating the drift length of photogenerated charges, we 
find that polaron recombination is not limiting the photocurrent for annealed devices. Additionally, we applied 
Monte Carlo simulations on blends of conjugated polymer chain donors with acceptor molecules in order to 
gain insight into the polaron pair dissociation. The dissociation yield turns out to be rather high, with only a 
weak field dependence. With this complementary view on dissociation and recombination, we stress the im- 
portance of accounting for polaron pair dissociation, polaron recombination as well as charge extraction when 
considering the loss mechanisms in organic solar cells. 
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The performance of polymer based solar cells has increased 
steadily over the years, culminating in 6% power conversion 
efficiency achieved recently ^ For a further optimization, a 
more fundamental understanding of the recombination pro- 
cesses in these devices is necessary, but a unified model has 
not been presented yet. Indeed,the dominant loss mechanism 
of the photocurrent remains a controversially debated issue. 

In order to understand the importance of loss mechanisms, 
we briefly describe the steps taken from absorbed light to 
flowing photocurrent, as shown in Fig. [T] For details, we re- 
fer to the book by Brabec et al.^ Light transmitted through the 
transparent anode of the solar cell, typically indium tin oxide 
on glass, is absorbed mainly in one constituent of the pho- 
toactive organic semiconductor blend, the conjugated poly- 
mer. Singlet excitons are generated. These neutral excita- 
tions can move by diffusion, but due to their high binding 
energy cannot be separated spontaneously at room tempera- 
ture. If they do not reach a donor-acceptor interface within 
their lifetime, they will recombine radiatively (a), sending out 
photoluminescence. Reaching the heterointerface, however, 
an electron transfer to the acceptor molecular yields a polaron 
pair with almost 100% yield (b), which can be separated — 
assisted by processes discussed later in the text — with a high 
yield. On the way of the separated charge carriers to their 
respective electrodes, charges can be trapped^, or bimolecu- 
lar recombination (c) can take place, although usually with a 
low rat^^. Also, the extraction of the photogenerated charges 
is influenced by the selectivity of the electrodes, namely the 
surface recombination velocity for electrons and holes^ ^. 

Coming back to the discussion on the dominant loss mech- 
anism, polaron pair dissociation or polaron recombination. In 
2004, Mihailetchi et al.^ described the experimental photocur- 
rent of a bulk heteroj unction solar cell completely by polaron 
pair dissociation, applying the Braun-Onsager model^ also 
considering a minor contribution of diffusion The field- 
dependence of the photocurrent was described relative to the 
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FIG. 1 : The steps from light absorption to the flow of photocurrent 
in a polymer solar cell, and some loss mechanisms: (a) Photolumi- 
nescence due to excitons not reaching the dono-acceptor interface 
during their lifetime, (b) Polaron pair dissociation or geminate re- 
combination, respectively, (c) Bimolecular polaron recombination. 



voltage Vq, where the photocurrent becomes zero. As the 
properties of organic semiconductors under dark and illumi- 
nation are not identical, this is not a physically well-defined 
quantity. Last year, Ooi et al.^ presented investigations of the 
photocurrent of polymer based solar cells by a pulsed tech- 
nique. They found a point of optimum symmetry, Vpos, rel- 
ative to which the field dependence of the photocurrent is al- 
most identical in forward and reverse voltage direction. They 
assigned this voltage to the built-in voltage. Although not 
performing a quantitative fit to the photocurrent, they assign 
the shape of the field dependent photocurrent to the important 
process of charge collection, implying a field independent po- 
laron pair dissociation, which is in contrast to the results of 
Mihailetchi et al.^. In both publication J^, the influence of 
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polaron recombination is neglected. Indeed, recent measure- 
ments of charge carrier loss processes in bulk heteroj unction 
solar cells, applying the versatile Photo-CELIV technique, 
photo-induced charge extraction by linearly increasing volt- 
age^^, show that the recombination rate is very low^ ^ 
The charge carrier dynamics are found to follow the temper- 
ature dependence of the bimolecular Langevin theor>^^, but 
with a reduced prefactoi^ Other reports point out that the 
order of the decay is larger than two, which i s expe cted for 
bimolecular processes, but between 2.5 a nd 3. Expla- 
nations for this behaviour are still debatecPI^^E^EIl, but the ex- 
perimental findings agree that the recombination rates found 
in bulk heteroj unction solar cells are low. Nevertheless, the 
question wether the photocurrent in bulk heteroj unctions can 
be described by one dominant process of polaron pair disso- 
ciation, polaron recombination, and polaron extraction, or if 
a combination of all three is necessary, has not yet been ad- 
dressed. 

In this paper, we will present experiments and simulations 
on polaron pair dissociation and polaron recombination, in or- 
der to contribute to this issue. We find that a unified model 
will have to consider dissociation, recombination as well as 
charge extraction simultaneously, although polaron recombi- 
nation can be mostly neglected for annealed state-of-the-art 
bulk heteroj unction solar cells. 

Polymer based solar cells where processed by spin coat- 
ing poly(3-hexylthiophene-2,5-diyl) (P3HT):[6,6]-phenyl- 
Cei butyric acid methyl ester (PCBM) blends made from a so- 
lution of (20mg P3HT+20mg PCBM) per ml chlorobenzene, 
on poly (3 ,4-ethylenedioxythiophene) :poly (styrenesulfonate) 
covered indium tin oxide/glass substrates. The active layer 
was about 105nm thick. Aluminum anodes were evaporated 
thermally. The annealed samples were subsequently treated 
for 10 minutes at 140°C. P3HT was purchased from Rieke 
Metals, PCBM from Solenne. All materials were used with- 
out further purification. The preparation steps were done in a 
nitrogen glovebox with attached thermal evaporation cham- 
ber. 

The samples were characterized in a He contact gas cryostat 
by current-voltage measurements as well as the photo-CELIV 
method. After a nanosecond laser pulse, the solar cell being 
held at the built-in voltage, the photogenerated charge carriers 
are extracted after a given delay time by an extraction voltage. 
The details of this method are described in Ref.^. The inset of 
Fig. [2] shows the schematic measurement cycle. The sample 
is held at the built-in potential. A nanosecond laser pulse is 
absorbed in the sample, generating excitons, and from these 
polaron pairs are formed due to the highly efficient charge 
transfer. During a variable delay time, the charge carriers can 
recombine, reducing the carrier concentration within the de- 
vice. After the delay time, two consecutive triangular volt- 
age extraction pulses — in reverse bias direction — are applied 
in superposition of the built-in potential, separated by a short 
waiting time. The first voltage pulse extracts all remaining 
photogenerated charge carriers, plus charges from doping and 
injection. As the latter two mechanisms are not in the focus 
of our investigations, the second voltage pulse helps to dis- 
criminate them from the photogenerated charges, which have 




0.0 0.5 1.0 1.5 2.0 

t [x10"%] 

FIG. 2: Photo-CELIV measurement of a P3HT:PCBM solar cell at 
260K with variable delay between laser pulse and charge extraction. 
For longer delay times, loss mechanisms reduce the signal magni- 
tude. The right extraction signal was recorded without laser pulse, it 
shows the contribution of injected charge carriers. The inset shows 
the schematic Photo-CELIV experiment: laser pulse, and-after a de- 
lay time td — two voltage pulses in reverse bias direction for charge 
extraction. In the experiment, the extraction pulse length was 800 /xs, 
and the waiting time between the two consecutive extraction pulses 
200 /is. 

already been drawn from the sample during the first pulse. 
By changing the delay time between laser pulse and charge 
extraction voltage pulses, a time resolution in analogy to a 
pump-probe experiment can be achieved, yielding the car- 
rier concentration vs. time dependence. Another outstanding 
feature of photo-CELIV is the simultaneous determination of 
charge carrier mobility and concentration of the extracted car- 
riers. Its almost only drawback is that, in contrast to transient 
photoconductivity^^, electrons and holes cannot be discrimi- 
nated. 

Retracing the derivation of the photo-CELIV analysis as 
published by Juska et al.^^ in 2000 due to the authors of 
Ref.l^, and comparing it to the original result for the mobility 
fi (Ref."^, Eqn. (5)), we found a somewhat different outcome. 
Solving the Riccatti equation numerically, we find 



^ 3AtLx(l + 0.21Aj7i(0))' 

with the original factor 0.36 replaced by 0.21. L is the device 
thickness, A is slope of the extraction pulse, and Aj/jo the 
height of the current extration peak relative to the dielectric 
contribution. In the normalised parameter space At^^^ vs. 
Aj/jo, the error of Eqn. ([T]) is below 2% for Aj/jo < 6, 
whereas the original evaluation is only in the range of < 70%. 
Thus, we use Eqn. ([T]) for the evaluation of all our CELIV 
measurements. 

An exemplary photo-CELIV measurement with variable 
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delay time between laser pulse and charge extraction for an 
annealed solar cell at 260K is shown in Fig. [2] The evalua- 
tion of the extracted charges over delay time yields the time- 
dependent carrier concentration n{t), which we could fit with 
the carrier continuity equation, 
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G-R = G- Cinp. 



(2) 



Here, the spatial derivative of the current is expected to be 
zero at the built-in potential. G is the generation term; with 
the laser pulse ait = 0, it is zero thereafter. We have shown 
here only a bimolecular recombination rate R, also assum- 
ing the electron and hole concentrations n resp. p to be equal. 
7 = g'/e • /i is the Langevin recombination parameter, with the 
elementary charge q, the effective dielectric constant e, and 
the sum of electron and hole mobility /i. C is a correction 
factor: We found that polaron losses in the solar cells can be 
described only by a reduced Langevin recombination. This 
works well for pristine devices; for annealed samples, a car- 
rier concentration dependent prefactor leads to an even better 
agreement, yielding a third order decay. The latter is probably 
related to delayed recombination due to trapping^^, however, 
it makes a difference only for long decay times not that rele- 
vant to steady-state solar cell operation, we concentrate on the 
reduced bimolecular recombination here. We found ( to be 
between 0.1 and 10~^, the details of which are described in 
Ref.^. At temperatures above 260K, the evaluation was made 
difficult by charge carrier injection, which can already be seen 
in Fig. [2] for the right extraction peak, which was recorded 
without laser excitation. 

In order to consider how relevant polaron recombination 
is for a working solar cell, we evaluate the mobility-lifetime 
product /ir. It is a measure of how efficiently charges can be 
extracted before recombination. It is used to calculate the drift 
length dc, which should exceed the device length in order to 
yield an efficient charge extraction. The drift length is given 
as 



prF, 



(3) 



where F is the electric field acting on the charges. In order to 
calculate an effective lifetime, the recombination rate has to 
be reshaped to R = n/r, 



(4) 



We point out that this effective lifetime depends on the car- 
rier concentration. Now, the mobility-lifetime product can be 
easily derived. 



(5) 



From the experimental photo-CELIV data, n, /i and ( can 
directly be extracted. Thus, using the experimental data, the 
temperature dependent /ir product can be calculated accord- 
ing to Eqn. ([5|. The data for temperatures above 260K could 
not be used due to dark injection currents disturbing the data 
analysis. Generally, the carrier concentration — needed for the 
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FIG. 3: The /ir product for the pristine and annealed bulk hetero- 
j unction solar cells. 



calculation of the effective lifetime — was taken at the mini- 
mum delay time of 5/iS. For the annealed sample, the carrier 
concentration has its maximum at these short times. This is 
also true for the pristine sample at 1 80K, but for higher tem- 
peratures, the carrier concentration is already decreased some- 
what at the minimum delay time. The resulting underesti- 
mated concentration (factor 2 at 260K) yields a corresponding 
overestimation of r. However, as the laser used in the photo- 
CELIV experiment generally leads to clearly higher light in- 
tensities as compared to solar radiation, a solar cell under op- 
erating conditions will have lower steady state carrier concen- 
trations, and the resulting lifetimes could even be longer. In 
Fig. [3j the /ir product for the pristine and annealed sample is 
shown. 

In order to estimate the collection distance, the extraction 
field F = {Vbi — V) / L still needs to be determined. L is 
the sample thickness, V the applied voltage. The compen- 
sation voltage of photo-CELIV measurements, at which the 
transient photocurrent becomes zero, corresponds to the built- 
in potential Vbi. For the pristine sample we determine a com- 
pensation voltage of 0.95V, for the annealed sample of 0.85V. 
From the current-voltage characteristics under illumination 
of lOOmW/cm^ (Fig.|4]), the corresponding open circuit volt- 
ages can be determined. At room temperature, they are 0.74V 
and 0.62V, respectively, and higher at lower temperatures as 
shown in Fig. [5] Therefore, we estimate an internal voltage 
Vbi — V of approx. 0.2V at open circuit for both samples, and 
0.95 resp. 0.85V at short circuit for the pristine and the an- 
nealed sample. The resulting collection distance, calculated 
using Eqn. ([3]), is shown in Fig.|6] Generally, dc is larger than 
the device thickness L = 105nm at temperatures above 200K, 
except for the pristine cell at open circuit conditions. This im- 
portant result is in agreement with our published macroscopic 
simulation of the influence of the charge carrier mobility on 
the solar cell performance^^, and also with the carrier concen- 
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FIG. 4: Current- voltage characteristics of (a) pristine and (b) an- 
nealed P3HT:PCBM solar cell in dependence on temperature. 



tration evaluated in Ref."^. Thus, we find that in state-of-the-art 
polymer solar cells, the nongeminate, bimolecular recombi- 
nation of electrons and holes is not limiting the performance 
under most conditions. 

In order to achieve a complementary view on the limiting 
factors in polymer solar cells, we also consider polaron pair 
dissociation. This was done by performing Monte Carlo sim- 
ulations of hopping transport in a gaussian density of states 
within a cubic lattice of 100 x 30 x 30 sites (with one hop- 
ping site per cubic nm). Some sites are considered electron, 
some hole transporting, in order to simulate a blend system. 
The width of the gaussian density of states distribution of 
the polymer donor was set to 75meV, and to 60meV for the 
fullerene acceptor, according to experimental results of ours. 
The dielectric constant of the blend was chosen to be 3.7. 
An electric field was applied along the long axis, it was as- 
sumed to be constant. For the other two directions, periodic 
boundary conditions were introduced. Coulomb interaction 
of the charge carriers as well as geminate and nongeminate 
recombination were accounted for. We extended the model 
in order to take intrachain transport along conjugated seg- 
ments of the polymer chains into account. Therefore, we in- 
troduced an effective conjugation length CL, consisting of 4 
or 10 monomer units, along which the charge transport — due 
to delocalization — is assumed to be instantaneous. The hop- 
ping process between the conjugated segments was calculated 
using the Miller- Abrahams jump rate. The results were av- 
eraged over at least 200 runs, with ten polaron pairs per run. 
Details of this approach are described elsewhere^^ . 

Assuming a perfect exciton dissociation, we just considered 
the charge transfer state: polaron pairs. They were generated 
on adjacent sites, the hole on a donor site and the electron on 
an acceptor site. During the Monte Carlo steps of hopping, 
they can either separate — assisted by the energetic disorder 
and the electric field — or recombine. The latter process hap- 
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FIG. 5: Open circuit voltage Voc (top), short circuit current density 
jsc (middle) and fill factor (bottom) of a pristine (dashed line) and 
annealed (solid line) P3HT:PCBM solar cell in dependence on tem- 
perature. The pristine samples levels off at low temperatures due to 
double diode behaviour, the annealed device is limited by the built- 
in potential. The corresponding current-voltage characteristics are 
shown in Fig.|4] 



pens with an effective recombination rate /ceff = rj^ , corre- 
sponding to the inverse effective lifetime of the polaron pair. 
A polaron pair was considered separated once one of its con- 
stituents reached an electrode. 

The dissociation yield, describing the successful separation 
of a polaron pair, is shown in dependence on the conjuga- 
tion length in Fig. |7] We considered two electric fields, 
and lO^V/m. The latter corresponds to the short circuit case 
in a lOOnm thick device with a built-in voltage of IV. Thus, 
the two fields approximately span the fourth quadrant of the 
current-voltage characteristics of a bulk heteroj unction solar 
cell under illumination. For both effective lifetimes consid- 
ered, 10ns and 10/iS, the field dependence of the dissociation 
process is not negligible, but weak. A dissociation yield be- 
ing comparable to the external quantum efficiency — having 
its maximum around 80% — can thus only be achieved when 
considering the effective conjugation length of the polymer 
chains. The fast intrachain process leads to an efficient charge 
pair separation process. This significant outcome implies that 
in a state-of-the-art polymer solar cell under working condi- 
tions, the field dependence of the polaron pair dissociation 
will reduce the fill factor, but only slightly. 

In accordance with our Monte Carlo simulations, we point 
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FIG. 6: Collection distance of pristine (dashed line) and annealed 
(solid line) P3HT:PCBM solar cell in dependence of temperature. 
The internal field under short circuit and open circuit conditions was 
used in Eqn. ([3]). The dotted line indicates the thickness of the bulk 
heteroj unction solar cell of 105nm. 
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FIG. 7: Monte Carlo simulation of the dissociation yield of a polaron 
pair in dependence of the conjugation length CL at T = 300 K. Two 
different fields are shown, OV/m and 10^ V/m, corresponding to the 
flat band case and the short circuit condition, respectively. 



solar cells is only due to bimolecular polaron recombination. 
The authors found that the same bimolecular recombination 
current dominates the dark current and the illuminated solar 
cell under open circuit. We note that their methods only con- 
sider polaron recombination currents, not polaron pair recom- 
bination. 

From our investigations, we find the photocurrent in or- 
ganic solar cells to be influenced by polaron pair dissociation, 
polaron recombination and also charge extraction. However, 
we see the need to distinguish between pristine and annealed 
P3HT:PCBM solar cells in discussing the limiting factors. In 
annealed samples, the performance is not limited by polaron 
recombination. We expect, however, a slightly field dependent 
polaron pair dissociation, which also influences the fill factor 
and thus the photocurrent. In contrast, pristine devices show a 
more severly limited photocurrent and fill factor as compared 
to their annealed couterparts. From our experiments, we see 
that the drift length particularly at low fields is not sufficient 
for a high yield of charge extraction. Thus, the photocurrent 
of pristine P3HT:PCBM solar cells is limited by both, polaron 
pair dissociation and polaron recombination, but in annealed 
devices mostly the influence of polaron pair dissociation is 
seen. The process of charge extraction needs to be accounted 
for in both, pristine and annealed devices. 

In conclusion, we investigated the polaron pair dissociation 
and polaron recombination in pristine and annealed bulk het- 
erojunction solar cells based on P3HT:PCBM blends by us- 
ing photo-induced charge extraction experiments as well as 
Monte Carlo simulations. From our experiments, we find that 
the collection distance is larger than the device thickness un- 
der most conditions, except for the pristine cell at open cir- 
cuit. Therefore, the photocurrent of state-of-the art polymer 
solar cells is not limited by the recombination of free po- 
larons. Concerning polaron pair dissociation, we found very 
high separation yields due to the delocalization of charge car- 
riers along conjugated segments of the donor polymer chains. 
This is accompanied by only a weak field dependence of the 
yield in the range relevant for solar cells under working condi- 
tions, with a significant amount of dissociation already at zero 
field. The latter is driven by the energetic and spatial disorder, 
but only made possible by delocalization along the extended 
polymer chains. Finally, the photocurrent of polymer solar 
cells can only be described accurately when polaron pair dis- 
sociation, polaron recombination as well as charge extraction 
are taken into account. This can only be done with macro- 
scopic simulations applying the appropriate physical models. 
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